Occurrence of Banned and Currently Used Herbicides, in Groundwater of Northern Greece: A Human Health Risk Assessment Approach

The presence of pesticide residues in groundwater, many years after their phase out in European Union verifies that the persistence in aquifer is much higher than in other environmental compartments. Currently used and banned pesticides were monitored in Northern Greece aquifers and a human health risk assessment was conducted. The target compounds were the herbicides metolachlor (MET), terbuthylazine (TER), atrazine (ATR) and its metabolites deisopropylatrazine (DIA), deethylatrazine (DEA) and hydroxyatrazine (HA). Eleven sampling sites were selected to have representatives of different types of wells. Pesticides were extracted by solid-phase extraction and analyzed by liquid chromatography. MET was detected in 100% of water samples followed by ATR (96.4%), DEA and HA (88.6%), DIA (78.2%) and TER (67.5%). ATR, DIA, DEA, HA, MET and TER mean concentrations detected were 0.18, 0.29, 0.14, 0.09, 0.16 and 0.15 μg/L, respectively. Obtained results were compared with historical data from previous monitoring studies and temporal trends were assessed. Preferential flow was the major factor facilitating pesticide leaching within the month of herbicide application. Moreover, apparent age of groundwater and the reduced pesticide dissipation rates on aquifers resulted of long-term detection of legacy pesticides. Although atrazine had been banned more than 18 years ago, it was detected frequently and their concentrations in some cases were over the maximum permissible limit. Furthermore, human health risk assessment of pesticides was calculated for two different age groups though drinking water consumption. In all examined wells, the sum of the HQ values were lower than the unity. As a result, the analyzed drinking water wells are considered safe according to the acute risk assessment process. However, the presence of atrazine residues causes concerns related with chronic toxicity, since ATR R values were greater than the parametric one of 1 × 10−6 advised by USEPA, for both age groups.


Introduction
Safe drinking water from surface-and ground-water is essential for human health, quality of life and socio-economic development of humanity and is a prerequisite factor for the human population [1]. The largest body of freshwater in the European Union quality standards exceedances in numerous European countries such as Spain, Slovenia, Hungary, Portugal and Germany, over last decade [24][25][26][27][28]. High concentrations of TER were detected, too. For example, the maximum detected concertation in groundwaters in the Águeda River Basin from Spain and Portugal was 2.5 µg/L [24]. Surprisingly, a serious lack of human health risk assessment studies in groundwater is observed in Europe. Therefore, the aim of this study was to investigate the groundwater quality, the presence and the persistence of these compounds (i.e., ATR and its metabolites, MET and TER), and to characterize their temporal and spatial variability in the studied aquifer, 18 years after atrazine's ban. Moreover, a chronic risk assessment of side effects on human health by consumption of contaminated drinking water was conducted. According to our knowledge, this study is the first attempt to correlate the human health risk assessment with the determination of pesticides residues in groundwater in Balkans and by extension in Greece.

Study Area Description
Sampling area was chosen taking into account the persistence of target pesticides in groundwater in Northern Evros, Thrace District, North Greece, as reported by previous studies and the extensive cultivation of wheat (2313.63 ha), cotton (2286.54 ha), sunflowers (1817.71 ha), maize (1695.18 ha) and beets (322.39 ha). Therefore, samples were taken from 11 sampling points, at Ardas Valley (third biggest arable area in Greece), divided in three Groups A sampling network of shallow groundwater was established during a previous study 20 years ago consisted of 4 experimental boreholes (Group A) [2]. Furthermore, 5 drinking water wells were included, which supply Orestiada town and local villages with potable water (Group B). In addition, two active irrigation wells were chosen (Group C) to include all available well types. Therefore, the sampling groups form was: Group A: A1, A2, A3 and A4, Group B: B1, B2, B3, B4 and B5 and Group C: C1 and C2. Sampling sites were located close to villages Rizia (A1, B1, C2), Keramos (A2), Plati (A3), Fylakio (A4, B2, C1), Elia (B3), Arzos (B4) and Kastanies (B5) (Figure 1). Slovenia, Hungary, Portugal and Germany, over last decade [24][25][26][27][28]. High concentrations of TER were detected, too. For example, the maximum detected concertation in groundwaters in the Águeda River Basin from Spain and Portugal was 2.5 μg/L [24]. Surprisingly, a serious lack of human health risk assessment studies in groundwater is observed in Europe. Therefore, the aim of this study was to investigate the groundwater quality, the presence and the persistence of these compounds (i.e., ATR and its metabolites, MET and TER), and to characterize their temporal and spatial variability in the studied aquifer, 18 years after atrazine's ban. Moreover, a chronic risk assessment of side effects on human health by consumption of contaminated drinking water was conducted. According to our knowledge, this study is the first attempt to correlate the human health risk assessment with the determination of pesticides residues in groundwater in Balkans and by extension in Greece.

Study Area Description
Sampling area was chosen taking into account the persistence of target pesticides in groundwater in Northern Evros, Thrace District, North Greece, as reported by previous studies and the extensive cultivation of wheat (2313.63 ha), cotton (2286.54 ha), sunflowers (1817.71 ha), maize (1695.18 ha) and beets (322.39 ha). Therefore, samples were taken from 11 sampling points, at Ardas Valley (third biggest arable area in Greece), divided in three Groups A sampling network of shallow groundwater was established during a previous study 20 years ago consisted of 4 experimental boreholes (Group A) [2]. Furthermore, 5 drinking water wells were included, which supply Orestiada town and local villages with potable water (Group B). In addition, two active irrigation wells were chosen (Group C) to include all available well types. Therefore, the sampling groups form was: Group A: A1, A2, A3 and A4, Group B: B1, B2, B3, B4 and B5 and Group C: C1 and C2. Sampling sites were located close to villages Rizia (A1, B1, C2), Keramos (A2), Plati (A3), Fylakio (A4, B2, C1), Elia (B3), Arzos (B4) and Kastanies (B5) (Figure 1). Figure 1. Study area of a transboundary aquifer (among Greece, Turkey and Bulgaria) and sampling sites. The capital letters A, B and C represent the experimental boreholes, drinking water wells and irrigation wells, respectively. Numbering identifies the villages close to sampling sites. Figure 1. Study area of a transboundary aquifer (among Greece, Turkey and Bulgaria) and sampling sites. The capital letters A, B and C represent the experimental boreholes, drinking water wells and irrigation wells, respectively. Numbering identifies the villages close to sampling sites.
Sampling was performed on a monthly basis for the period 28 March-5 September 2018, conducting 5 samplings. The samples (volume of 1 L) in triplicate were poured in dark polypropylene bottles. Then, they were transported to the Laboratory of Agricultural Pharmacology and Ecotoxicology, Department of Agricultural Development, at Democritus University of Thrace in ice-boxes and were stored in the dark at temperatures below −20 • C until analysis that was performed within a period no longer than 1 week. Experimental borehole samples were manually pumped using an experimental tube. The drinking water and irrigation wells were equipped with a pump system and samples were collected automatically before chlorination stage.

Chemicals and Instrumental Analysis
The pesticide standards had the highest available purity (>97%) and were purchased by Dr. Ehrestrofer GmbH (Augsburg, Germany). Physicochemical properties of studied compounds are shown in Table 1. The HPLC grades, acetonitrile, ethyl acetate, water and methanol for liquid chromatography were purchased by Riedel de Haen (Seelze, Germany). LiChrolut ® EN Polymer-based solid-phase extraction cartridges with 200 mg absorbent and 3 mL volume were purchased by Merck (Darmastdt, Germany). Groundwater samples were prepared for instrumental analysis using solid phase extraction (SPE) for the multi-residue analysis as described by Papadakis and Papadopoulou-Mourkidou [30] with slight modifications. Water samples of 1 L were extracted by cartridges which were preconditioned with adding of 4 mL methanol followed by 4 mL deionized water. Samples were passed through cartridges at a flow rate of 5 mL/min. Target compounds were eluted with 7 mL methanol followed by 3 mL ethyl acetate. Next, samples were concentrated under nitrogen stream at 50 • C. Finally, samples were dissolved with 1.25 mL of methanol and stored at −20 • C until instrumental analysis.
Samples were analyzed by a HPLC/DAD equipped with autosampler (Finnigan Surveyor, Thermo Scientific, Boston, MA, USA). The analytical column C18 Speedcore 100 × 4.6 mm was purchased by Fortis Technologies Ltd. (Cheshire, UK). Chromatographic data were processed by the ChromQuest 5.0 software (Finnigan Surveyor, Thermo Scientific). The mobile phase was consisted of acetonitrile (A) and water (W). The flow was set at 1.0 mL/min and the gradient included the following steps: the elution began at 20-80/A-W, 20-80/A-W (0-20 min), 95-5/W-A (20-25 min), 95-5/A-W (25-26 min) and 20-80/A-W (26-33 min). Total run time was 40 min. The injection volume was 25 µL. The column oven temperature was adjusted at 30 • C. MET, TER, DEA and DIA were detected at 220 nm, while ATR and HA at 240 nm. For further confirmation of the target peaks, the UV absorption spectra taken at the apex of each sample were compared with those obtained from the standard solutions and control spiked samples. The quantification was performed using external working standard calibration curves (1, 10, 50 and 100 µg/mL). The accuracy (recovery) and precision (repeatability) of the analytical method were evaluated with the analysis of fortified (at 0.1 µg/g and 0.5 µg/g) tap water samples in duplicate. The limits of detection (LOD, µg/L) were determined as the lowest concentrations giving a response of three times the baseline noise of the analysis of three control samples. The limits of quantification (LOQ, µg/L) were determined as the lowest concentrations of a given compound in fortified samples that could be quantified with relative standard deviation lower than 20%. Positive detections of ATR, DIA, DEA, MET and TER were also confirmed with gas chromatographic analysis using a Trace 2000 gas chromatograph connected with the GCQ plus ion-trap mass spectrometer (Thermoquest, Austin, TX, USA). Gas chromatographic analysis was carried out on a 30 m × 0.25 mm I.D., 0.25 µm film thickness CP-SIL 8 CB (5% phenyl, 95% dimethylpolysiloxane) low bleed/MS column (Varian Analytical Instruments, Middelburg, The Netherlands) and the GC and MS operational conditions were those mentioned by Vryzas et al. [31]. Fortified tap water samples were used for the validation of the analytical method. Fortification levels includes the LODs, LOQs and the 0.1 µg/L. For all compound the LODs were ranged from 0.001 to 0.005 µg/L and LOQs from 0.01 to 0.05 µg/L. The recoveries were higher than 86% for all compounds with relative standard deviation lower than 15% at tested fortification levels.

Human Health Risk Assessment
Human health risk assessment was conducted for atrazine, metolachlor and terbuthylazine. According to Kim et al. [32] human health risk assessment of pesticides can provide information about the probability and the kind of effects to human population. In this case, oral exposure through drinking water consumption was considered as pathway to people. Risk assessment was divided to carcinogenic and non-carcinogenic and to two age groups, adults and children. Drinking water was provided to the local population by wells located close to the villages Elia, Arzos, Fylakio, Rizia and Kastanies ( Figure 1).

Chronic Daily Intake (CDI)
CDI shows the estimated intake amount of pesticide per kilogram body weight (Equation (1)): where D IP is the average daily intake, EF i is the exposure frequency (365 days per year for both age groups), ED i is the exposure duration (6 and 70 years for children and adults, respectively), BW i is equal to 70 kg for adults and 20 kg for children and AT is the average lifespan (2190 and 25,550 days for children and adults, respectively). The determination of the average daily intake (D IP ) was estimated using the Equation (2), which was suggested by Muhammad et al. [33], Papadakis et al. [4] and Ali et al. [34]: where C i (µg/L) represents extreme and mean concentration of pesticide residues and IR i shows the intake rate of water (0.87 L/day for children and 1.41 L/day for adults).

Hazard Quotient (Non-Carcinogenic Risk Assessment)
To calculate the hazard quotient (HQ), CDI was divided with the respective reference dose of each compound (Equation (3)): where RfD is the acute toxicity reference dose [35].
The RfD values for atrazine, metolachlor and terbuthylazine, were 0.035, 0.015 and 0.008 (mg/kg-day), respectively [36,37]. When HQ values are equal or greater than 1, the exposed part of population is under health risk.
Multiple pesticides residues risk (HQs) can be calculated by the sum of HQ for individual pesticide using the Equation (4):

Carcinogenic Risk Assessment
Carcinogenic risk (R) was calculated by the Equation (5) [4,32]: where SF is the cancer slope factor (mg/kg-day), which reflects the possibility of the individual pesticide to cause cancer and ADAF is an age factor considering the early life pesticide exposure (3 for children and 1 for adults). Among the studied pesticides the only available SF value is for atrazine with value 0.22, provided by IRIS [36].

Statistical Analysis
An analysis of variance was conducted (ANOVA) for human health risk results using a mixed linear model according to Mas et al. [38]. The calculations were performed using the Excel Solver Add-in Package.

Concentrations and Detection Frequency
The sampling sites were in the Ardas valley, an aquifer vulnerable to pesticide contamination according to previous studies [5,23]. In general, the mean concentrations presented a wide range of values exceeding, occasionally, the maximum permissible limit of 0.1 µg/L in experimental boreholes (group A), drinking water wells (group B) and irrigation wells (group C) ( Table 2). According to the results shown in Table 3, the ATR mean concentration was above 0.1 µg/L in A2 and A4 sites. DIA mean concentration was found higher than 0.1 µg/L in A4, while HA surpassed it in A2. DEA mean concentration were higher than 0.1 µg/L in A1 and A4. TER presented limit exceedances in A2 and A4 and MET in all boreholes. As far as drinking water wells, ATR showed higher mean concentration than 0.1 µg/L in B1, B2 and B5. DIA mean concentration was higher in B1, B3 and B5. In HA case, mean concentration was above the limit in all drinking water wells. The mean concentration of MET was higher than 0.1 µg/L in B5, while TER limit exceedances were found in B1, B2 and B3. Concerning the irrigation wells, exceedances of 0.1 µg/L were observed for ATR, and DEA in both wells, for MET in C2 and for HA, DIA and TER in C1. Furthermore, the highest concentration of ATR, DIA, DEA, HA, MET and TER were 1.86, 2.91, 0.65, 0.46, 0.93 and 1.0 µg/L, respectively (Table 3). Table 3. Summary of positive samples in pesticide residues and exceedances of the EU permissible limit.  TER  ATR  DIA  DEA  HA  MET  TER   A1  5  80  80  60  60  80  60  0  40  20  20  40  0  A2  5  80  80  80  80  80  40  0  20  0  0  20  40  A3  5  80  80  40  80  80  20  40  20  0  0  40  20  A4  5  80  80  40  80  80  60  40  40  20  20  20  60  BI  5  80  100  100  80  100  100  60  60  40  0  40  100  B2  5  80  100  100  80  100  60  60  40  0  0  0  60  B3  5  100  100  100  100  100  80  0  100  40  40  20  0  B4  5  80  80  80  80  80  40  0  40  0  0  0  20  B5  5  100  100  100  100  100  80  100  100  80  40  60  80  C1  5  100  100  100  60  100  80  60  100  60  40  60  80  C2  5  100  100  100  100  100  80  60  40  60  40 Table 4 presents the percentage of positive samples in pesticide residues and the percentage of samples with concentration higher than permissible limits set by EU. Throughout the study, DIA and MET presented 90.9% detection frequency in water samples followed by ATR (87.3%), DEA and HA (81.8%) and TER (63.6%). The highest exceedance frequency of the maximum permissible limit in wells was observed in well B5 followed by C1, C2, B1 and B3. The measured values of pesticide concentrations in all sampling sites are presented by the box and whisker plot in Figure S1.    [41] found that DIA was the most frequent detected compound (72%), followed by TER (50%), with maximum concentrations 1.42 µg/L for DEA, 0.40 µg/L for DIA and 0.46 µg/L for TER. Moreover, a third Spanish study in agricultural areas showed maximum concentration 0.33 µg/L for ATR, 0.37 µg/L for DEA, 0.34 µg/L for TER and 0.55 µg/L for MET with detection frequency ranged from 4% (DEA) to 68% (MET). The same study presented results from Portuguese groundwater in agricultural areas, too. TER had the highest concentration (1.88 µg/L) with frequency detection reached 56%, followed by ATR (0.19 µg/L) and detection frequency 25%. DEA and MET concentration were lower than 0.1 µg/L [24]. According to Korosa et al. [27], in groundwater samples from Slovenia ATR and DEA were detected at concentrations up to 0.228µg/L and 0.10 µg/L and their frequency of detection was 94.6% and 98.2%, respectively. In another study conducted in United Kingdom and France, the highest concentrations in British groundwaters for ATR, DIA, and DEA were 0.20, 0.10 and 0.16 µg/L, respectively, while the highest concentrations were lower than 0.1 µg/L, in France [42]. Table 4 presents the respective data of similar monitoring studies.

Historical Vulnerability of the Transboundary Aquifer to Contamination by Pesticide Residues
The fact that target compounds reached concentrations above the quality standard values for drinking water indicates that studies made during pesticides registration process are not always complying with the results from monitoring studies. It is estimated that less than 1% of the pesticides applied reach the target pest and the remaining distributed to various environmental compartments including groundwater bodies [43]. In a previous study was found that waters of that site constitute an aquatic mixture with different residence times and various leaching mechanisms are involved to the pollution of groundwater [44].
Similarly to the results obtained 15 to 19 years ago extreme concentrations were observed in this study, indicating the presence of point source pollution sites nearby the studied wells. Target compounds had been monitored previously (between 1999-2003), at the same locations, before atrazine ban in EU [2]. Moreover, a similar study was conducted between 2010-2012 (data not shown), confirming the occurrence of ATR, DEA, DIA and MET (provided by the Greek Ministry of Agricultural Development). In order to have a better perspective on pollution temporal trends, the data of this study were compared with those of 1999-2003. Fifteen to nineteen years ago MET had been detected at least once in 63 % of the wells followed by ATR (61%), DEA (50%), alachlor (47%) and DIA (34%). According to Vryzas et al. [2], maximum concentrations for ATR (1.48 µg/L), DEA (0.76 µg/L), DIA (0.07 µg/L) and MET (1.54 µg/L) had been detected at the same drinking water wells sampled in this study and considerable higher pesticide concentrations were detected in shallow groundwater from experimental boreholes. This study confirms previously reported data on adsorption, degradation and leaching of selected pesticides [2,31,[44][45][46]. Vulnerability of the aquifer to pollution depend on the land uses, soil properties, geological characteristics of the unsaturated zone, the hydraulic properties, the depth of the vadose zone and the leaching potential or physicochemical properties of the contaminant.
According to previous studies focused on this area, ATR degraded faster than MET in all soils of the vadoze zone and the biotransformation rates of both compounds decreased as the soil depth increased. Hence, the chronic presence of ATR in field is indicated by the higher biotransformation rate of ATR in soil taken from the middle of a studied field in comparison with soil sampled from the field margins [2]. The major metabolites of ATR and MET were found at higher concentrations in the 10-20 cm layers of all soil cores studied 0-110 cm below ground surface. However, the enhanced biodegradation rate of ATR in these soils is not enough to prevent the contamination of groundwater bodies. Our results have been observed by other studies, too. According to McMahon et al. [47], Kolpin et al. [48] and Steele et al. [49] degradation rates of triazine parent compounds are slower than their transport rates in groundwater. Moreover, Vonberg et al. [25] and Jablonowski et al. [50] conducted a solid samples analysis contaminated with atrazine from a lysimeter and agricultural field, respectively, reporting the persistence of the parent compound for over 20 years after its last application, followed by possible leaching.
Adsorption studies of ATR, DEA, DIA, HA and MET were also conducted in soils from five different depths (i.e., 0-10, 10-20, 20-40, 40-80 and 90-110 cm below ground surface). The present study revealed that when pseudo-equilibrium stage reached, the amount of the adsorbed compounds accounted only for 10, 14, 27, 43 and 94% of the initial amount of DEA, DIA, ATR, MET and HA, respectively, spiked to the soils. According to this study, it was expected that more than 57 and 73% of the applied dose of MET and ATR, respectively, to be desorbed into the soil water and be available for leaching to deeper soil layers [45].
In addition, the low adsorption capacity of ATR and MET within soil profile of the studied area, proved that the preferential flow is a major pesticide leaching mechanism in this area since pollutants can reach the saturated zone of the aquifer through preferential flow paths (shrinkage of the clay minerals, plant roots and earthworms forming burrows) without going through chromatographic flow within the unsaturated zone and thereby circumventing the degradation processes [46]. Studies on the apparent age of the studied aquifers shown that the residence time of groundwater bodies ranged from 1.2 to 50 years [2].
The leaching mechanisms prevailed in this area has been also studied in an extensive four-year field experiment focused on soil water samples taken from 0-25, 35, 60, 100 and 160 cm below ground surface [46]. According to this study, MET, ATR, DEA and DIA were detected in more than 67% of the total soil water samples. The main conclusion of this study was that the corn-applied herbicides have been leached below the surface soil via macropore-dominated pathways in less than one month after their application. Agricultural practices (application of pesticides and sprinkler irrigation) used in this area, soil structure and hydrogeological conditions increase the leaching potential of pesticides in the studied area. It is worth notice that alachlor, another banned herbicide, with very limited half-life period (DT 50field = 14 days) had been detected in soil water of the studied area at concentrations greater than 0.1 mg/L up to 40 months after its application.
Irrigation is usually carried out by self-propelled sprinkler irrigation systems or basin irrigation systems multiple times in summer period. These types of irrigation system provide high volumes of high-pressure water, which, in combination with rainfall, can exacerbate the phenomenon of leaching [46]. Moreover, Nouma et al. [51] mentioned that basin irrigation systems can favor pesticides leaching to groundwaters, too. As recommend by Vryzas, et al. [46], the late pesticide application, use of drip instead of sprinkler irrigation and delayed first irrigation seem to be the major management actions according to good agricultural practice that prevent pesticide leaching to groundwater in a semiarid Mediterranean region. The limited spatial and temporal variation of concentration levels observed in studied wells indicates a continuous load of the aquifer with the target compounds. The continuous use of high amounts of atrazine for more than 30 years was enough to contaminate the soil and aquifer and to be detected with its metabolites in groundwater 15 years after its last use (2004). However, illegal applications cannot be excluded since the studied area is 20 km from Turkish and Bulgarian borders and illegal trade of banned pesticides had been observed.
The variability and seasonality of temperature and rainfall impact significantly on the fate of pesticides within soil zone. Along the study, the temperatures were ranged between 5.3 and 39.9 • C, reaching extremely high values. The vapor pressures of pesticides can be increased as the temperature increases. In this study, pesticides have low vapor pressures ( Table 2) and therefore volatilization could not cause losses, except for DEA. Moreover, high rainfall values were observed in May and June, 48.6 and 59 mm, respectively ( Figure 2). Heavy rainfall could induce the remobilization of organic pollutants adsorbed in soil or sediments, especially lipophilic compounds such as ATR, TER and MET. Moreover, high rainfall intensities can cause overflow events and increase the frequency of leaching phenomenon, which can also be more rapid under rainy conditions [52,53]. The MET and ATR effectiveness against corn weeds and the limited available herbicides resulted their extensive use during the period of 1980-2005. ATR withdrawal brings out the TER as the most used herbicide, until nowadays. Application rates are similar, ranged from 600 to 1200, 900 to 1800 and 860 to 1200 gr/ha for ATR, MET and TER, respectively. The cropping system and major crops has been gradually changed from 2005 till now. However, field crops are still the major crops in the area and the irrigation practices are the same that were used 20-40 years ago. TER applied by farmers in agricultural area of Ardas valley as a pre-emergence herbicide for maize, cotton and beet cultivations in The MET and ATR effectiveness against corn weeds and the limited available herbicides resulted their extensive use during the period of 1980-2005. ATR withdrawal brings out the TER as the most used herbicide, until nowadays. Application rates are similar, ranged from 600 to 1200, 900 to 1800 and 860 to 1200 gr/ha for ATR, MET and TER, respectively. The cropping system and major crops has been gradually changed from 2005 till now. However, field crops are still the major crops in the area and the irrigation practices are the same that were used 20-40 years ago. TER applied by farmers in agricultural area of Ardas valley as a pre-emergence herbicide for maize, cotton and beet cultivations in April, replacing the banned ATR. The extended treatment period of TER coincides with the significant rainfall period in the studied area. The use of a basic measure also identifies TER as a moderate leachers: in fact, TER has a GUS (groundwater ubiquity score) score of 2.19, indicating that it is a leacher [54].
In laboratory investigations with various top-soils, the degradation of [ 14 C] TER was examined to determine which soil features would enhance metabolite production. After a 45-day incubation period, the concentration of TER dropped by 30%, while desethylterbuthylazine was generated at the same time. In different soils, the TER half-life ranged from 88 to 116 days. In addition, a research team hypothesized that TER is more readily absorbed by soil than atrazine, which could shield the molecule against biodegradation [55].
Commercial plant protection products bentazone, terbuthylazine and metolachlor were used to investigate the degradation and leaching of bentazone, TER and S-metolachlor, as well as their metabolites applied to a weighable, monolithic, high precision lysimeter with a loamy, sandy soil. Bentazone decomposed faster in soil than TER and S-metolachlor, but TER metabolization after the second application resulted in higher formation of desethylterbuthylazine and TER hydroxylation. On day 12, the greatest concentration of TER (310 g/kg) was found while the topsoil contained 16 g/kg of terbuthylazine on day 150, indicating leaching losses [56].

DEA to ATR Ratio
The DEA to ATR ratio (DAR) has been used to categorize point and non-point source pollution of groundwater and, in order to, characterize the degradation and transport of ATR in response to its metabolite DEA. This ratio can give us an indication of the major leaching mechanisms contribute to the pollution of groundwater and the capacity of the unsaturated zone to biodegrade ATR to DEA. During the transport of atrazine through chromatographic flow within the biological more active unsaturated zone it could be metabolized in significant amounts by microorganisms to DEA [2,57,58]. In such cases the DAR would have values higher than 0.4, or even close to 1.0. In contrary, when atrazine bypasses the vadose and enters the saturated zone through preferential flow the contact time between atrazine and soil microbial community could be shorter and, therefore, the DAR ratio would be less than 0.4. The DAR ratio can provide information about atrazine leaching behavior based on the fact that atrazine represents a closer adsorption capacity to DEA than to HA in spite of HA was found as the main metabolite of atrazine at same area. In addition, this soil can adsorb higher amount of atrazine than DEA. Therefore, DEA can be leached faster than atrazine through chromatographic or preferential flow [46]. The calculated DAR in this study was found to be higher than one in some cases, and lower than one in most samples ( Table 5), indicating that contamination in some cases comes from diffuse sources but most probably the bound atrazine was gradually desorbed from the soil matrix to the soil water and moved to groundwater through preferential flow [25,46,59]. Although DAR values can be used to characterize the predominant leaching mechanisms in the studied aquifer, the different types of wells (experimental boreholes, drinking water wells and irrigation wells) used in this study and the respective differences in construction techniques restrict the estimation of the various mechanisms involved in the leaching processes. Our results are in agreement with those of Vryzas et al. [44], conducted in the same area 15-19 years ago, who found similar DAR values few months after the application of atrazine. Overall, atrazine's degradation products showed higher concentrations than did the parent compound. DIA exhibits a large range of concentrations varying between 0.01 µg/L and 2.91 µg/L. According to biotransformation studies conducted in the soil profile of the studied area HA was the most frequently detected metabolite and with the highest concentrations. The second most frequently detected degradation product in soil was DEA, followed by rare DIA detections [2]. The overwhelming majority of soil water samples with DEA presence, showed DEA had greater concentrations than DIA and the ratio values CDEA/CDIA reached 33 [51]. DEA (50% of groundwater samples) was more frequently detected than DIA (34% of groundwater samples) in an extensive groundwater monitoring program conducted in the same area 15-19 years ago [2]. Contrary to previous reported data, in this study, atrazine and its metabolites were detected with equal frequency of detection.

Human Health Risk Assessment
An extended discussion was preceded related to the presence, occurrence and distribution reasons of target pesticides at studied area. Results on human health risk assessment are presented in Table 6. Although, the HQ values for individual pesticide did not exceed the value 1, the estimated non-carcinogenic risk for children was higher, when compared to adults. The HQ values for mean pesticides concertation were ranged between 0.0171 to 0.1913 for adults and between 0.0393 to 0.5752 for children. The highest mean values are reported to metolachlor and the lowest to atrazine. The highest HQ values were determined in drinking water well close to Rizia namely, 0.2507 and 0.7817 for adults and children, respectively. Similar HQ values for atrazine and metolachlor in surface water were reported by Papadakis et al. [4], in a Greek study. The risk level for TER was low with HQ values lower than 0.6.
The sum of HQ values did not reach the unity in all studied wells. The greatest cumulative potential risk was determined in the Rizia well with values 0.2836 and 0.8299 for adults and children, respectively. The lowest potential risk has the Elia well, with values lower than 0.4. Consequently, according to the acute risk assessment, the studied drinking water wells are characterized as safe.
Oppositely, the carcinogenic risk assessment showed high values. In all cases, ATR R values were higher than the parametric one of 1 × 10 −6 recommended by USEPA, for both age groups, showing that the local population were under carcinogenic risk ( Table 6). The water consumption through the Fylakio well presents the highest risk, while the Arzos well was the lowest. The R values are ranged between 0.0002-0.0018 for adults and 0.0012-0.0181 for children. Previous study conducted in Greece, indicated high carcinogenic risk only for children [4].

Conclusions
All pesticides were detected in both shallow and deep ground-water bodies (experimental boreholes, drinking or irrigation water wells). DIA and MET presented the highest detection frequency in water samples, while TER and DIA showed the highest exceedance frequency of the EU permissible limit. Perennially, these herbicides were used in a variety weed-sensitive cultivations, including wheat, beet, cotton, maize and sunflower. The type of crops grown in a region largely determines the types of pesticides that were used, which in turn affects the residues of pesticides that are typically found in the groundwaters of studied region. Moreover, illegal used of banned pesticides cannot be excluded from a transboundary area. Despite agricultural use of atrazine has been banned in Greece for more than 18 years ago, atrazine and its metabolites residues are still detected in groundwater of the region, indicating their high persistence in saturated zone. The repeated application of studied pesticides could lead to enhanced biodegradation, as previously reported in the studied area, but the remaining amounts of bound residues was gradually desorbed from the soil matrix to the soil water and moved to groundwater through preferential or chromatographic flow. Furthermore, the detection of extreme concentrations suggests the presence of occasional point-sources pollution. However diffuse sources cannot be excluded. These observations were confirmed by DAR ratios and the frequency of detection of parent compounds and their respective metabolites, suggesting multi-mixed leaching. The drinking water consumption for local people is safe considering the acute risk assessment. However, the atrazine R values suggested high carcinogenic risk. We suppose that EU regulatory agencies could benefit from the linking of monitoring data with probabilistic human models in order to propose efficient pollution management methods.  Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest. The article reflects only the author's view and the Agency is not responsible for any use that may be made of the information it contains.